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Abstract. High statistical accuracy of experiments KAS- 
CADE and ANI allowed to obtain approximations of primary 
energy spectra and elemental composition in the "knee" re- 
gion. Obtained results point out to the correctness of QGS JET 
interaction model and 2-component model of primary cosmic 
ray origin up to 100 PeV energies. 



1 Introduction 

Absolute differential EAS size spectra around the knee mea- 
sured at different atmosphere depths and different zenith an- 
gles are not explained yet from the point of view of a single 
A — AAir interaction model and a single model of primary 
energy spectra and elemental composition. Such an attempt 
has been made in work (Ter-Antonyan and Haroyan, 2000) 
based on an unified analysis of KASCADE, AKENO, EAS- 
TOP and ANI EAS size spectra. The results of approxima- 
tions of primary energy spectra by rigidity-dependent steep- 
ening spectra pointed to the correctness of QGSJET interac- 
tion model and two-component composition of primary pro- 
ton spectrum in the knee region. 

Here, on the basis of KASCADE (Glasstetter et al., 1999) 
and ANI (Chilingarian et al., 1999) EAS size specti-a the 
multi-component model of primary cosmic ray origin (Bier- 
mann, 1993) has been tested in the framework of method 
(Ter-Antonyan and Haroyan, 2000). 



2 Method 

The testing of primary energy spectra was carried out in 0.03- 
500 PeV primary energy range using -minimization (Ter- 
Antonyan and Haroyan, 2000) 
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is detectable EAS size spectra measured in i = 1 , ... to size 
intervals and j = 1, . . .n zenith angular intervals (see Fig.l, 
symbols), is an energy threshold of detected EAS elec- 
trons, N*^{E > i?e) is the estimation value of EAS sizes 
obtained by detected electron lateral distribution functions at 
observation level (t); 
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are the expected EAS size spectra at d'^A/dEo energy spec- 
tra of primary nuclei (A = 1, . . . , 59); 
^ = m-n— p— lisa degree of freedom at p number of 
unknown parameters. 

The function We in expression (3) is determined in general 
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where dQ./dNf. is an EAS size spectrum at the observation 
level [t) for given Eo,A,9 parameters of a primary nucleus 
and depends on A — AAir interaction model; 

Ae — cos 01 — cos 02; 
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is a probability to detect an EAS by scintillation array at EAS 
core coordinates \x\ < X/2, \y\ < y/2 and to obtain estima- 
tions of EAS parameters (N*, s - shower age, x*,y* - shower 
core location) with given accuracies; 

d^/dN* is a distribution of N*(Ne, s, x, y) for given EAS 
size (iVe). 
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In most of EAS experiments the EAS cores are selected 
in P ~ 1 range providing a log-Gaussian form of a measur- 
ing error {d'^ JdN*) with an average value \n{Ne ■ 6) and 
RMSD (In, where 5 involves all transfer factors (an energy 
threshold of detected EAS electrons, 7 and ji contributions) 
and slightly depends on Eq and A. In these cases, one may 
standardize the measured EAS size spectra to the EAS size 
spectra at the observation level 



dl(0,e,t) dI(Ee,0,t) 



dNe 



dm 



(6) 



where r/ = exp{(7e - l)'^a%/2} and 7e is the EAS 

size power index. 

Taking the above into account, the expected EAS size spectra 
(3) can be estimated according to the expression 
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where the sum is performed into limited number (fcmax) of 
nuclear group. 

2.1 2-component primary energy spectra 

Energy spectra of primary nuclei = 1 — 59) according 

to the multi-component model of primary cosmic ray origin 
(Biermann, 1993) are presented in 2-component form: 
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where the first component (ISM) is derived from the explo- 
sions of normal supernova into an interstellar medium with 
expected rigidity-dependent power law spectra (Biermann, 
1993) 
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and the second component (S W) is a result of the explosions 
of stars into their former stellar winds with expected rigidity- 
dependent power law spectra (Biermann, 1993) 
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where ^a is a scale factor {E in TeV units) from approxima- 
tions (Wiebel-Sooth and Biermann, 1998); 



EiSM = RiSM ■ Z, Es 
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are the corresponding rigidity-dependent cut-off energies of 
ISM-component and knee energies of SW-component; 
RiSM and Rsw are model parameters of magnetic rigidities 
of corresponding components and Z is the charge of A nu- 
cleus. 

The values of model predictions (Biermann, 1993) for spec- 
tral parameters are: 

71 = 2.75 ± 0.04, 72 = 2.67 ± 0.03, 73 = 3.07 ± 0.1 

and rigidities Rism ^ 120 TV, Rsw - 700 TV at factors 
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Fig. 1. KASCADE EAS size spectra at different zenith angles 
(Glasstetter et al., 1999) (symbols). Lines correspond to expected 
EAS size spectra according to QGSJET interaction model and 2- 
component origin of primary cosmic rays. 



of uncertainty ~ 2. 

The fractions of each component iSA,i = Si{EA,A), i = 
1, 2) are determined according to 



Sa,i = 1 - Sa,2 

Sa,2 = {2ZRism)^^'-^^''^ 



(12) 
(13) 



at 7A=i,o = 2.75 and 7^>i,o = 2.66 (Biermaim, 1993). The 
expressions (12,13) are consequences of normalization of (8- 
10) to approximation of balloon and satellite data (Wiebel- 
Sooth and Biermann, 1998) at = 1 TeV. 

Thus, minimizing -functional (1) on the basis of mea- 
sured values of dI{6i)/dNej and corresponding expected 
EAS size spectra (7) at given m zenith angular intervals and 
n EAS size intervals one may evaluate parameters of primary 
spectra for given (fcmax) nuclear groups. Evidently, the accu- 
racies of solutions for spectral parameters strongly depend 
on the number of measured intervals (m • n), statistical errors 
and correctness of dQ{Eo , A, 9, t) / dN^ determination in the 
framework of a given interaction model. Moreover, the value 
of points out a reliability of applying primary model. 
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Fig. 2. Expected primary all particle and P, He, O, Fe energy 
spectra obtained by approximation of KASCADE data. The 
symbols are CASA-BLANCA (Fowler et al., 2000) and AGASA 
(Yoshida et al., 1995) data. 



Differential EAS size spectra dfl{Eo, A, 0, t)/dNe for given 
Eo = 0.032, 0.1,..., 100 PeV, A = 1,4,12,16,28,56, t = 
0.5, 0.6, . . . , 1 Kg/cm^, cos 9 = 0.8, 0.9, 1 were calculated 
using CORSIKA562(NKG) EAS simulation code (Heck et 
al.,1998) at QGSJET (Kalmykov and Ostapchenko, 1993) in- 
teraction model. Intermediate values are calculated using 4- 
dimensional log-linear interpolations. Estimations of errors 
of expected EAS size spectra dfl/dNg at fixed Eo,A,d,t 
parameters did not exceed 3 — 5%. 
Unknown parameters in minimization (1) were: 
Eis M - cut-off energy of ISM proton component; 
Esw - knee energy of proton SW-component; 
73 - power index after the knee (Esw)'- 
77 - a systematic shift (discrepancy) from expression (6). 
Moreover, power indices 71,2 were changed too in the range 
of the model uncertainty and the relative uncertainties of ex- 
pected spectra (ap) in the -minimization (1) we set equal 
to 3%. 



3 Results 

The testing of primary model predictions by minimizations 
(1) were carried out on the basis of KASCADE (Glasstetler 
et al., 1999) (t =1020 g/cm^) and AMI (ChiUngarian et al., 
1999) (700 g/cm^) EAS size spectra at 5 zenith angular inter- 
vals. The KASCADE EAS size spectra (symbols) and cor- 
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Fig. 3. Expected fractions of SW-component at different nuclei ob- 
tained by approximation of KASCADE EAS size spectra. 



responding expected spectra by 2-component model (lines) 
are shown in Fig.l. The obtained primary all-particle energy 
spectrum and energy spectra of P, He, O, Fe nuclear groups 
are presented in Fig. 2 in comparison with CASA-BLANCA 
(Fowler et al., 2000) and AGASA (Yoshida et al., 1995) mea- 
surements. The obtained values of approximation parameters 
of primary energy spectra are: 

71 = 2.78 ± 0.03; 72 = 2.65 ± 0.03; 73 = 3.28 ± .07; 
EisM = 210 ± 60 TeV; Esw = 1900 ± 100 TeV; 
and J] = 1.03 ± 0.03 at = 1-22. 

Obtained fractions of the primary stellar wind (SW) compo- 
nent 



AswiEA,A) = 6A^ 



d^2/dEA 
' d^/dEA 



(14) 



versus energy at different primary nuclei are presented in Fig. 
3. The results are extrapolated up to 0.1 TeV primary energy 
range. 

The testing of 2-component model of CR origin was also car- 
ried out by ANI EAS size spectra (Chilingarian et al., 1999) 
measured at the mountain level. The results are shown in 
Fig. 4 and correspond to x^ = 1-05 and systematic shift 
T} = 1.14. The systematic underestimations of KASCADE 
(77 = 1.03) and ANI (77 = 1.14) EAS size spectra (see expres- 
sion (6)) one may explain by energy thresholds of detected 
EAS electrons {E^ ~ 3 MeV for KASCADE and Eg ~ 10 
MeV for ANI experiments). 

The values of spectral parameters obtained by approxima- 
tions of ANI data agree with corresponding parameters ob- 
tained from approximations of KASCADE data except for 
cut-off energy Eism = 460 ±100 TeV at ANI data analysis. 
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Fig. 4. ANI EAS size spectra at different zenith angles (Chilingar- 
ian et al., 1999) (symbols). Lines correspond to expected EAS size 
spectra according to QGSJET interaction model and 2-component 
origin of primary cosmic rays. 



4 Conclusion 

Predictions of 2-component model of the cosmic ray origin 
and QGSJET high energy interaction model allow us to ex- 
plain the measured EAS size spectra in the knee region with 
accuracy better than 10%. 

All model parameters of 2-component primary energy spec- 



tra, obtained from EAS data agree with theoretical predic- 
tions (Biermann, 1993) in the frames of standard errors with 
the exception of 73 parameters. Obtained spectral slopes 
(73 = 3.28 ± 0.07) after the knee of SW-component are sig- 
nificantly steeper than 2-component model predictions ((Bier- 
mann, 1993), 3.07 ± 0.1) and this result requires of further 
investigations. 
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